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ABSTRACT 

The shear triplet statistics is a geometric method to measure cosmological parameters with 
observations in the weak gravitational lensing regime towards massive haloes. Here, this pro- 
posal is considered to probe the dark energy equation of state and its time derivative in view 
of future wide-field galaxy surveys. A survey with a median redshift of ^ 0.7 and a total 
area of ^ 10000 square degrees would be pretty effective in determining the dark matter 
cosmological density and in putting useful constraints on the dark energy properties. 
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1 INTRODUCTION 

The deep theoretical understanding of gravitational lensing makes 
it an attractive probe of dark energy, one of the main puzzle of mod- 
ern cosmology. Dark energy can show up either in a pure geometri- 
cal way by affecting the distance redshift relation or via its effect on 
the growth of structure. Lensing is sensitive to both signatures since 
a ratio of distances appears in the scaling of lensing parameters 
with redshifts and because the effective mass of lensing structures 
reflects the po wer spectrum and growth rate of large-scale density 
perturbations jSartelmann & Schneide3l200ll ; iMunshi et ai]|200d 
and references therein). Probing dark energy through both geome- 
try and growth by taking weak gravitational lensing two-point func- 
tions (such as the shear power spectrum) of d istant galaxy images 
as a function of redshift can be appealing HHeavens et al J 120061 : 
lHuteredl2003 : ISong & Kno?3l2004l : iTakada & Jainll2004 and ref- 
erences therein), but this method relies on the interpretation of the 
distortion signal on scales where non linear evolution necessitates 
accurate modeling from numerical simulations. In the strong lens- 
ing regime the distance ratio and the cosmological parameters can 
be derived using observations of giant luminous arcs in rich clus- 
ters of galaxies but the mass profile and the mass normalization 
of the lenses must be assumed to be known, for examp le from X- 
ray measurements ( ISerenol2002l : ISereno & Longdl2004l , and refer- 
ences therein). Strong gravitational lensing observations of image 
separations an d time delays could provide useful constraints as well 
( lLindeill2004t) . 

The basic idea behind geometric methods is to separate out 
the information purely from the distance ratios, irrespective of the 
lensing mass distribution. In the strong lensing regime, this can be 
achieved by considering multiple images systems and the ratios of 
arc angular positions, but this still requires the assumption that th e 
lensing mass, even if unknown, is pretty regular dOolse et alj2002h. 
The weak lensing regime seems more promising. Ijain & Tavlod 



( l2003l) proposed to take the ratio of the galaxy-shear correlation 
function at different redshifts behind galaxy groups and galaxy 
clusters. Dark energy parameter constraints based on this geometric 
method, which extends into the nonlinear matter power spectrum 
but still drops out the mass of the lens, have been then discussed by 
several authors in view of future galaxy sur veys ( Bernstein & Jaig 
|2004| : IHu & Jair]|2004l : ISong & Rnoxll2004l) . The method was then 
further developed by directly considering ratios of the shears be- 
hind a cluster without fi rst generating the cross-correlati on func- 
tions dTavlor et al.l2007 ^. Alternativelv. lZhang et ai] ( l2005h pointed 
out how the cross-correlation of a foreground galaxy-density field 
or shear field with the shear field from a background source pop- 
ulation scales with the source redshift in a way that can be used 
to constrain cosmology without making assumptions about the 
mass/galaxy power spectrum. 

In this paper, we want to reconsider the so called triplet statis- 
tics, an original idea to constrain cosmolog ical paramete r s from 
weak lensing in galaxy clusters proposed in iGautret et all ( l2000t) . 
This method is able to disentangle the effect of the lensing mass, 
described by local convergence and shear terms, from the cos- 
mological parameters by considering the ellipticities in triplets of 
galaxies located at about the same angular position bu t having dif- 
feren t redshifts. Differently from similar proposals dTavlor et al.l 
|2007|) . the triplet statistics is not limited to the out skirts of mas- 
sive ha loes and can be used even in the inner regions. lOautret et al.l 
( l2000h originally considered the method for determining the cos- 
mological constant with observations towards very massive galaxy 
clusters. Here, we review the triplet statistics and discuss the mea- 
surement of the dark energy properties with future weak lensing 
survey. 



2 BASICS 
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The distortion of images of background galaxies is determined 
by the convergence fc, i.e. the lensing strength, and the complex 
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shear 7 = 71 + 172- The lensing parameters can be related to 
the value they would have for a source at a reference redshift 
dSeitz & Schneideilll997l : rGautret et alfcOOOh . 



k = Uj{Zs)krct, 

7 = aj(2;s)7,rf. 



(1) 
(2) 



The lensing factor uo(zb) is defined as a ratio of angular diameter 
distances and contains the cosmological dependency. 



(3) 



where Dij is the angular diameter distance between the redshifts 
Zi and Zj with the redshifts of interest being those of the ob- 
server o, the lens d and the source s. We consider a standard 
Friedmann-Lemaitre-Robertson-Walker model of universe, filled 
with non-interacting pressure-less matter and dark energy, param- 
eterized through its equation of state w ~ Wo + Wa{^ — a), with 
a = l/(l + z). Since the contribution from relativistic particles is 
negligible in the redshift range investigated in our analysis, we will 
neglect it in what follows. In such a model of universe, the angular 
diameter distance between an observer at Zi and a source at 2, is 



Hol + zj IOkoI^/^ 
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dz 
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- nva-3(i+"'o+»a)e-3«^a(i-a) ^ nKa-2, (5) 



where 

Hjz) 
Ho 

and Ho is the present value of the Hubble parameter; Qm and fiv 
are the today normalized densities of dust and dark energy, respec- 
tively; r^K = 1 — fiivi — f^v; Sinn is defined as being sinh when 
Ok > 0, sin when JIk < 0, and as the identity when JIk ~ 0. For 
the exp ression of the distance in an inhomogeneous universe, we 
refer to lSereno et al] j200lll2002l) . 

The transformation from intrinsic to observed ellipticity in 
the weak lensing regime (7, fc ^ 1) tak es a simple form 
( ISeitz&SchneideJll997l ; lOautret et alj|2000l) . Due to a lensing 
hal o, a galaxy with int rinsic ellipticity Cs is imaged with elliptic- 
ity iG^treTrtiDllQoO) 



e - (1 -ff^)es +g 



(6) 



with g = 7/(1 — fc) being the complex reduced shear. Background 
galaxies at (nearly) the same angular position probe the same local 
cluster mass distribution and gravitational potential. Then compar- 
ing the shear amplitude for three galaxies having different redshift 
allows to separate the effect of the mass distribution from cosmol- 
ogy. A geometrical operator can be built from the measured ellip- 
ticities and redshifts Zi of the three galaxies i — \a, b, c| in such 
a way that it depends only on cosmology CGautret et aU200Q) 



Tabc — 



1 LOa iOaeb^c 
1 OJi, LObecCl 
1 LJc UJceaCl 



(7) 



Redshifts of the galaxies inside triplets are in ascending order. For 
three intrinsically spherical galaxies, Tate ~ when the lensing 
factors LUi are computed for the actual values of the cosmological 
parameters, apart from noise. Tate is linear regarding to the ellip- 
ticities, which makes the principal source of noise randomly dis- 
tributed around zero. Apart from noise, the main part of the triplet 
operator contains the cosmological dependence and can be approx- 
imated as 
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Figure 1. Confidence regions in the Cm — f^A plane as expected for a 
lensing survey with Zm ~ 0.7. Contours show the 1,2-3(t limits around the 
fiducial model Qm = 0.3, Qa = 0.7. The inner thin and the outer thick 
contours refer to a survey with total area A = 10000 dcg^ and 1500 deg^, 
respectively. 
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where the apex denotes that the angular diameter distances have 
been calculated for the actual values of the cosmological param- 
eters and where we have considered only one component of the 
shear. If we neglect the contribution from the local convergence, the 
reduced sher g can be identified with the local shear and then it is 
enough to take the ratio of the observed ellipticities of a pair of near 
galax ies in order to separate the effect of cosmology dTavlor et al.l 
I2OO7I) . Neglecting k introduces a systematic error that can be signif- 
icant for the largest clusters. A 10% variation of the cosmological 
parameters changes the reduced shear by about 1% which is ten 
times small er than the relative variation due to the 1 — tjfcrof for 
fcrcf ^ 0.1 dOautret et alj|2000l) . Even if the majority of the signal 
comes from intermediate mass clust ers, the highest ret urn in accu- 
racy comes from the largest haloes jTavlor etal ] |2007l) so that this 
effect must be properly accounted for. 

Let us sort the redshift in a triplet such that Za < z^ < Zc- 
From the matrix form of the operator Tabc, it is clear that if two 
galaxies in a triplet are very close, then Tabc — with no regard 
to the cosmological parameters. Tabc — also if the minimum 
redshift in the triplet is very close to the lens redshift. Once fixed 
the minimum and the maximum redshift in a triplet, the sensitiv- 
ity of the operator Tabc is maximized for an intermediate redshift 
Zb nearly in the middle of the redshift range. Once fixed Za and 
Zb, the sensitivity increases with the maximum redshift Zc- Hence, 
the main information from the triplet method comes from the high 
redshift tail of the background source distribution. 
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Figure 2. Confidence regions in the — plane for a lensing survey 
covering 10000 deg^. Contours show the 1,2-3(T Umits around the fiducial 
fiat model Qm = 0, ioq = —1 and assuming no evolution for the dark 
energy, Wa = 0. The dashed fine separate models with either accelerated 
or decelerated expansion. 



3 FORECAST FOR LENSING SURVEYS 

There are a number of current and planned imaging surveys for 
weak lensing analyses jPeacock et alj llOO^ . Beyond 2007, the 
funded VST ( VLT Survey telescope) public survey KIDS will cover 
at least 1500 square degrees in 4 broadband optical filters. Com- 
bined with a following near infrared coverage by VISTA, this will 
yield a 9-band optical-IR survey with depth approximately 2 mag 
deeper than the Sloan and with very accurate photometric redshift 
estimate. The typical KIDS lenses will be at z ~ 0.2 — 0.4. Be- 
yond KIDS, the planned Dark Energy Survey on the CTIO Blanco 
telescope or the darkCAM survey at VISTA are expected to take 
another steps forwards in terms of sky coverage, imaging ~ 10000 
square degrees. In what follows, we give a cosmological parameter 
estimation forecast for such surveys. 

The background redshift distribution for a typical magnitude- 
limited survey can be taken to be 



dz 



= no 



3£^ 

2zl 



exp 



(9) 



The redshift scale zq is related to the median redshift of the sur- 
vey, Zm, by z(i — Zin/1A12; no gives the total number density of 
sources with usable photometric redshift and shape estimate. For a 
KIDS-like survey we can take Zm ~ 0.7. 

The main part of the lensing signal comes from halo lens- 
ing masses in the range 5 x 10^^ < M/Mq < 10^^ The num- 
ber density of haloes can be accurately calculated and varies from 
^hai ^ 10^ per square degree for M ^ 10^^ to riha i ^ 10~^ 
per square degree for M ~ 10^^ jTavlor et al. I [2^. An aper- 
ture size dmax — 2' — 5' corresponds roughly to the virial radii of 
such massive haloes over the redshift range considered. An inner 
circle with aperture 9in ~ 0.1' — 0.8' has to be excised from the 
data to exclude the arclets and the strong lensing regime. Consid- 



Figure 3. Confidence regions in the wq — Wa plane for a lensing survey 
covering 10000 dcg^. Contours show the l,2-3o- limits around a reference 
flat model wg = —1, uiq = assuming flat geometry and a sharp prior 
Qm = 0.3. 



ering a simple isothermal sphere model for the lenses, we get that 
(7rof) ~ 0.04 — 0.15 in the mass range of interest. We remark 
as the triplet operator is proportional to {■y^^{) rather than {7rcf)^- 
These estimates for the shear are lower limits since considering a 
Navarro-Frank-White profile as a deflector, the shear signal would 
increases for the same halo mass. 

In order to study the sensitivity of the triplet method on the 
cosmological parameters, we perform a analysis. The main dif- 
ference with the statistical analysis in Gautret et al. (2000) is that 
they considered a mean triplet operator as the average over all 
triplets whereas here we consider a built on the linearly inde- 
pendent triplets. Not all of the triplets we can put together behind 
a lens contain independent information. It can be easily shown that 
from Ng near galaxies only Ng — 2 triplets out of the binomial fac- 
tor {Ng , 3) are independent. If a triplet contains at least one galaxy 
not already included in the sample of other triplets, then it is lin- 
early independent. Triplet selection can be properly optimized. As 
maximum distance between the triplet components we take a sep- 
aration of Ar ~ 20" (Gautret et al. 2000). This is a good balance 
between having a sufficient number of triplets and not smearing the 
signal. As x^, we consider 



X 



E 

i ,{a ,b,c} 



/ ^abc 
V 5Tabc 



(10) 



with the sum running over the lensing haloes I optically selected in 
the survey and the independent triplets for each halo. Foreground 
lenses and background sources are modeled according to the pre- 
vious discussion. For each background galaxy we completed the 
triplet by selecting the two neighboring galaxies (Ar < 20") 
which maximize the signal. 

Statistical and syste matics errors affecting the method have 
been deeply discussed in iGautret et all ( [2000i) . The main sources 
of noise are: i) the intrinsic source ellipticities; ii) the errors on 
measured ellipticities; in) the fact that sources do not experience 
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exactly the same potential and finally iv) the errors on measured 
(photometric) redshifts. It can be shown that due to the linearity 
of the operator, the noise is linear with respect to each individual 
term and then is proportional to 1 /\/iV with N the total number of 
triplets. The dominant contribution to the error budget comes from 
the intrinsic ellipticity. 

Together with statistical noise, several systematics affect the 
method. The main ones are well understood and have been identi- 
fied as i) a bias due to an asymmetry in the probability distribution 
of the terms Aa; due to photometric redshift errors and, mainly, ii) 
contamination by background structure, either galaxy-galaxy lens- 
ing or large-scale structure. Sources in very close angular pairs, for 
which another galaxy could play the role of lens, could be rejected 
from the sample. It can be shown that for surveys large enough the 
effect of large-scale structur e can be neglected at first order with re - 
spect to the statistical noise dOautret et al.l200(]| : lTavlor et al.l2007h . 
The dominant noise sour ce is by due to intr insic ellipticity disper- 
sion dOautret et alj|200d iTavlor et al.ll2007h . Then, as an estimate 
of 5Tabc in Eq. dlOt we will take 

where the factor 2 in the denominator arises because we are us- 
ing only one component of the measured ellipticity. As intrinsic 
ellipticity dispersion we take CTe = 0.3. The effect of lensing by 
large-scale structure will be discussed later However, even if we 
are underestimating the statistical noise by considering only the 
main contributor, the choice of parameters has been generally con- 
servative. Results for a survey with median redshift ~ 0.7 and a 
typical density of no ~ 30 galaxies per square arcminute are shown 
in Figs.[T][2]and[3] As a first step we have considered dark energy 
in the form of a cosmological constant {wo = — 1 and Wa = 0), 
see Fig.[T] The contours show the 1, 2 and 3cr confidence limits for 
two parameters (Ax^ = 2.30, 6.16 and 11.8 respectively). With a 
survey area of 10000 square degrees, the triplet statistics can con- 
strain the matter density parameter but is pretty insensitive to the 
total amount of vacuum energy. This degeneracy changes with red- 
shift. For tests in the strong lensing regime based on luminous gi- 
ant arcs, the contours are nearly orthogonal to the triplet estimator 
since the distance ratio is evaluated for very high redshift sources 
dSereno & Longoll2004h . 

Assuming a flat model and a constant equation of state, the 
triplet statistics can put a quite firm upper limit on wo, see Fig. [2] 
On the other hand, the confidence regions spreads out well within 
the phantom regime {wo < —1), as usual for methods based on 
the distance-redshift relation. The phantom regime would be also 
consistent with an higher value of the matter density parameter. 
Contours in the wq — Wa plane are still pretty elongated, see Fig. [3] 
On its own, the triplet statistics can not say much about the evo- 
lution with time of the dark energy, but if combined with orthog- 
onal methods such as the cosmic microwave background radiation 
or the baryonic acoustic oscillation in the matter power spectrum, 
constraints could be significant. 

The linear size of the confidence regions shrinks approxi- 
mately as 1/ Ng, with Ng the total number of galaxies in the 
survey. Cosmological constraints then would strongly benefit by 
a very large survey area. Furthermore, increasing the median red- 
shift of the survey would both increase the local galaxy density and 
would probe the distance ratio in a redshift range more sensitive to 
cosmology. 

Together with the variance term proportional to the intrinsic 



uncertainty per shear mode due to the galaxy intrinsic ellipticities 
{(7^), which is related to the shape noise and to the shot noise, the 
other main source of error is due to lensing by large-scale struc- 
ture in between the len s and the sources ((tlls). related t o the 
sampling variance term dZhang et alj[20o3 : lHu & Jainll2004l) . De- 
pen ding on the survey strategy, these two terms can be compara- 
ble. iTavlor et a i] l l2007t) gave a simple approximate scaling rela- 
tion between the two terms for a survey of median redshift Zm, 
a'lss — {24:.lz^Az)a'^, where Az is the typical photometric red- 
shift error for the survey. Then, for 2;m ~ 0.7 and Az ~ 0.05, error 
estimates on cosmological parameters obtained considering just the 
shot noise should be increased by ~ 15%. 



4 FINAL REMARKS 

The aim of this paper has been the evaluation of the triplet statis- 
tics as a dark energy probe in view of future galaxy survey. The 
main source of statistical uncertainty in our statistical approach was 
the intrinsic source ellipticity. As regards the main systematics, the 
contribution of large-scale structure to the observed shear, being 
uncorrelated with the cluster effect, should av erage out over inde - 
pendent clusters along different lines of sight dTavlor et al.ll2007h . 
In a survey measuring photometric redshifts, data should be col- 
lected in redshift bins with width equal to the typical redshift error 
at that redshift. Then, the number of independent triplets behind a 
cluster would be Nun — 2, with Nbin the total number of bins. 
However, photometric redshifts, especially if some infrared filters 
are available, should be accurately determined. 

In this paper, we have considered any halo mass profile but the 
method could be optimized using the symmetrical properties in the 
mass distribution of galaxy groups and clusters. A nearly elliptical 
matter distribution would allow to consider tangential shear aver- 
aged in concentric annuli, i.e. to collect triplets selecting galaxies 
in the concentric ring instead of a small local patch. This would 
make nearly sure that for any galaxy we can find a pair of galaxy 
redshifts that maximize the signal. 

Some conclusions on the viability of the shear triplet method 
can be drawn by comparison with the shear ratio geometric test. 
The confidence regions we plotted seem larger that those ob- 
tained with the shear ratio test in a survey with similar properties 
(Tavloretal. 2007). As we have seen the shear test is biased for 
large mass haloes M > 10^^ M©, where the reduced shear should 
be properly considered instead of the shear. Being these haloes 
pretty rare if compared with halo masses of order of < IO^^Mq, 
which provide the bulk of the signal, this systematic effect should 
not jeopardize the shear ratio method. In any case, since the two 
techniques require the same kind of measurements they should be 
properly integrated. This is desirable especially because the analy- 
sis of optically selected large mass halo in wide field survey should 
begin with very massive haloes which, on turn, pay the highest div- 
idend. 
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